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^ ■ ABSTRACT 

\) , Context. In about half of Seyfert galaxies, the X-ray emission is absorbed by an optically thin, ionized medium, the so-called "Warm 

• . Absorber", whose origin and location is still a matter of debate, 
i-^ ' Aims. The aims of this paper is to put more constraints on the warm absorber by studying its variability. 

Mh' Methods. We analyzed the X-ray spectra of a Seyfert 1 galaxy, Mrk 704, which was observed twice, three years apart, by XMM- 

f~^ , Newton. 

^ ■ Results. The spectra were well fitted with a two zones absorber, possibly covering only partially the source. The parameters of the 

absorbing matter - column density, ionization state, covering factor - changed significantly between the two observations. Possible 
J^ ' explanations for the more ionized absorber are a torus wind (the source is a polar scattering one) or, in the partial covering scenario, 

an accretion disk wind. The less ionized absorber may be composed of orbiting clouds in the surroundings of the nucleus, similarly 

to what already found in other sources, most notably NGC 1365. 

►-^ . Key words. Galaxies: active - X-rays: galaxies - X-rays: individual: Mrk 704 

in 
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^^ 1. Introduction The paper is organized as follows: in Section 2 we report 

'""' I on the XMM-A^evvfo« observations and data reduction, while 

• , Absorption from ionized matter in the X-ray spectrum of AGN ^j^^ ^^^^^^^ ^^^^ ^^^^ -^ ^^ discussed in Section 3. Section 4 
^ (the so-called warm absorber) was discovered many years ago ^^^^^^^ ^^^ ^^^^ -^ ^^ ^^^ ^ ^^-^^ ^^^ ^^ ^^ observations, while 
^ (Halpern 1984). Since then, many advances in its understanding ^^^ ^^^^^^^ ^^ summarized and discussed in Section 5. 

,__! have been made, and we know now that it is present in about 
r* half of Seyfert galaxies (e.g. Reynolds 1997) and that the mat- 

.^ ter is photoionized. Warm absorbers are also known to vary, 2. XMM-Newton: Observations and data reduction 
K> and indeed the first discovered absorber was variable (Halpern 

j_| ■ 1984). The location of the warm absorber (or, indeed, of the ab- XMM-A^ewfo« observed Mrk 704 twice, on 2005-10-21 (OBSID: 

C3 ■ sorbers, as more than one ionizing zone is often found) is, how- 0300240101) and on 2008-1 1-02 (OBSiD: 0502091601). In both 

■ " " ' ever, largely uncertain. There is some evidence for its origin as cases, EPIC pn and MOS were in Small Window mode (apart 

a wind from the dusty torus envisaged in unification models for from MOS2 in the first observation, which was in Full Frame 

Seyfert galaxies (Blustin et al. 2005), but cases in which an ori- mode), which ensures that no significant pile-up is present, as 

gin from the disk seems to be preferred do also exist (Krongold verified with the EPATPLOT tool, in the pn detectors (the only 

et al. 2007). oi^c used for the analysis). Mrk 704 is by far the brightest source 

Mrk 704 is a local (z=0.029234) Seyfert 1 .2 galaxy ( Veron- i" the field of view. 
Cetty & Veron 2010), bright enough in X-rays to be detected by Data were reduced with SAS 10.0.0, using caUbration files 
SwiftfBAT (Ajello et al. 2008). In this paper we report on ex- generated on 2010-6-1 1. Screening for intervals of flaring par- 
treme warm absorber variability, on yearly time scales, revealed tide background was done consistently with the choice of ex- 
by two XMM-A^evv'fon observations, and possible variability on traction radii, in an iterative process based on the procedure to 
monthly time scales from short Swift/XRT observations. maximize the signal-to-noise ratio described by Piconcelli et al. 
(2005). After this process, the net exposure time was of about 15 



Send offprint requests to: G. Matt, e-mail: and 68 ks for the 2005 and the 2008 observation, respectively, 
matt @f is .uniromaS . it adopting extraction radii of 36 and 40 arcsec, and patterns to 4. 
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Fig. 1. The 0.3-2 keV (upper panel), 2-10 keV (middle panel) and the hard-to-soft ratio (lower panel) light curves, for the first (left 
panels) and second (right panels) XMM-Newton observations. 



The background spectra were extracted from source-free circular 
regions with a radius of 50 arcsec. The same regions were used 
for the timing analysis. Spectra were binned in order to over- 
sample the instrumental resolution by at least a factor of 3 and 
to have no less than 30 counts in each background-subtracted 
spectral channel. The latter requirement allows us to use the x^ 
statistics as a goodness-of-fit-test. RGS source and background 
spectra were extracted with standard procedures, adopting the 
data reduction pipeline RGSPROC, and choosing the NED opti- 
cal nucleus of Mrk 704 as the reference point for the attitude 
solution. 



3. XMM-Newton: Data analysis 

We will use in this paper only data from the EPIC pn and the 
RGS camera, because after the latest release of the EPIC pn re- 
distribution and the RGS contamination model (May 2010) their 
cross-calibration is as good as 3% (A. Pollock, private commu- 
nication). 

In Fig. [U the 2-10 keV (upper panel), 0.3-2 keV, (middle 
panel) and the (2-10 keV)/(0.3-2 keV) hardness ratio (HR) light 
curves are shown for the first (left) and the second (right) obser- 
vations. Intra-observation variability in the count rates is appar- 
ent in both cases. While in the first observation no clear variation 
in the HR is present, in the second observation an anticorrela- 
tion of the hardness ratio with the flux is evident, indicating that 
the source gets softer when brighter. Therefore, while for the 
first observation we will use, for the spectral analysis, the time- 
integrated spectrum only, for the second observation we will also 
make a time-resolved spectral analysis. 

Spectral analysis has been performed using the XSPEC 
package, version 12.6.0. Unless specifically stated, all eiTors re- 
fer to the 90% confidence level for one interesting parameter. 

Let us start with the time-averaged analysis of both obser- 
vations. In Fig. |2]flie best fit models and residuals are shown, 
assuming a simple power law absorbed by the Galactic column 
density in the direction of the source (2.97x10^" cm^", Kalberla 
et al. 2005). The fits are clearly bad, not surprisingly given the 
extremely simple model. From the residuals, a soft excess, a 
broad absorption trough between about 1 and 2 keV - indica- 
tive of obscuration by ionised gas - and an iron line at 6.4 keV 
are apparent. 

Comparing the spectra of the two observations (Fig. [3j, a 
huge difference in the soft part of the spectrum is clear, suggest- 



ing a much stronger absorption in the first observation. Above 
a few keV, instead, the two spectra are pretty similar (see also 

Fig. nil. 

3.1. The iron line 

We first concentrate on the iron line, and to this end we limit 
the bandwidth to 5-10 keV, where the soft excess and the warm 
absorber do not have an important impact on the spectral model- 
ing. More specifically, the choice of 5 keV as the lower limit of 
the band is due, on one hand, on the necessity to have the entire 
line within the band (which is verified a posteriori), and on the 
other hand to avoid any effect of the warm absorber In fact, a fit 
with a simple power law gives very different values - especially 
for the first observation - of the photon index depending on the 
lower limit of the band, stabilizing only above 5 keV. The pho- 
ton index steepens increasing the lower limit, a clear indication 
of absorption at low energies. 

We fitted the 5-10 keV spectrum with a power law, a 
Compton reflection component (PEXRAV model in XSPEC, pa- 
rameter R - the solid angle subtended by the reflecting matter to 
the nuclear source in units of In - fixed to 1) and an iron line, 
firstly modeled with a naiTow Gaussian at 6.4 keV {cr fixed to 
0). Results are summarized in Table[Tl where best fit spectral pa- 
rameters for this and following models are reported. In both ob- 
servations, the statistical quality of the fit is not very good, and 
residuals are evident around the line, suggesting it is resolved. 
Indeed, leaving cr free to vary, the quality of the fit significantly 
improves (see Table [T]), with cr of 100-150 eV, corresponding to 
a FWHM velocity 2-3 times larger than that of the Broad Line 
Regions (about 5000 km/s, Stirpe 1990). Letting the centroid en- 
ergy of the line free to vary, instead, does not improve the quality 
of the fit. 

The most likely broadening mechanism is therefore mat- 
ter rotation in an accretion disk. We therefore substituded the 
Gaussian with the DISKLINE model (valid for Schwarschild met- 
ric) with outer radius fixed to 1000 gravitational radii and the 
emissivity index to -2.5. The rest-frame line energy is again fixed 
to 6.4 keV. The fits are of comparable statistical quality to those 
obtained with a broad Gaussian. The moderate width of the line 
translates in a quite large inner radius (much larger than the in- 
nermost stable orbit) and in a small inclination angle. 

On the other hand, it is possible that the line is actually a 
blend of a narrow line (originating e.g. in the torus) and of a 
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Fig. 2. The spectrum and residuals for the first (left) and second (Right) XMM-Newton observations when fitted with a simple 
power law absorbed by Galactic interstellar matter. 
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Fig. 3. The two XMM-Newton spectra shown together for com- 
parison. The soft X-ray flux in the first observation is much lower 
than in the second observation, while above 5 keV the two spec- 
tra are almost indistinguishable. 



relativistic line. We therefore fitted the spectrum with the com- 
bination of the two, fixing for simplicity the inner radius of the 
relativistic line to 6 gravitational radii (the innermost stable orbit 
for Schwarschild metric) and the inclination angle to 30 degrees. 
The fit is again good, and we cannot really distinguish between 
the different models on statistical ground. (A fit of similar qual- 
ity is obtained if a LAOR model is used for the line, again with 
the inner radius fixed to the innermost stable orbit.) We limit 
ourselves to note that a relativistic component seems to be re- 
quired, but whether the line emitting region is truncated or a nar- 
row component is also present cannot be told. 

Comparing the two observations, there is marginal evidence 
for a stronger (and broader) broad component of the line in the 
first observation (see Fig.|4]i. The quality of the data is, however, 
not good enough to further explore this possible line variability. 

For both observations, and for all line models, the power law 
photon index is around 1.7-1.8. 

3.2. The warm absorber 

We now study the broadband (0.5-lOkeV) spectra of the source. 
Having in mind the residuals to a simple power law model 



^^"!"""i^^^ 
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Fig. 4. The hard X-ray spectra of the two XM-Newton observa- 
tions shown together for comparison. 



(Fig. [2|l, we added to the model already used to study the iron 
line (adopting for simplicity a broad Gaussian for the line) a 
second power law to model the soft excess, and one or more 
warm absorber region. The warm absorber was modeled with the 
ZXIPCF model, based on the XSTAR photoionization code. The 
main limitations of the model are a fixed ionizing continuum (a 
power law with photon index of 2) and fixed element abundances 
(equal to solar). The ionization parameter is ^ = L/nR^. The re- 
sults are summarized in Table |2] For the first observation, the 
addition of an emission line around 0.9 keV (to be likely identi- 
fied with the Neon IX He-like triplet) was also necessary, as also 
suggested by the RGS analysis. In fact, in the RGS spectrum of 
obs 1 a line at 0.903+0.001 keV (corresponding to the forbidden 
line of the triplet) with a flux of 1.69(+0.72)xl0"^ ph cm"^ s"' 
was found (see below). 

3.2.1 . The first observation 

Let us first discuss the first observation, where the absorption 
is much stronger. A single ionization zone is insufficient, the 
addition of a second zone improving significantly, even if not 
dramatically, the statistical quality of the fit. Moreover, in the 
single zone model F, (the photon index of the soft power law) 
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Fig. 5. Obs 1: Best fit model and data/model ratio for the two zones model. Left panel: full covering absorbers. Right panel: partial 
covering absorbers. 
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Fig. 6. Obs 1 : Contour plots (ionization parameters vs covering factors) for the low (left panel) and high (right panel) ionization 
absorbers. The three curves represent the 67%, 90% and 99% confidence levels. 



is unrealistically steep, while F/, (the photon index of the hard 
power law) is much flatter than the value found in the hard band. 

The addition of the second zone improves the quality of the 
fit, which however remains not fully satisfactory. Visual inspec- 
tion of the data/model ratio (see Fig. |5]l, however, shows that 
this is due mostly to noise rather than to unfitted features, even if 
some wiggles are present around the Neon line. Indeed, an im- 
provement (Ay^=- 16) of the fit is obtained by letting free to vary 
the width of the line, when a value of 48 e V is found. This may 
suggest a blend with other lines (mostly iron L), too faint to be 
individually detected by the RGS. However, it is also possible 
that the broadening of the line is an artifact to correct for a too 
strong, in the model, iron UTA (this may occur e.g. if the iron 
abundance is lower than solar). Unfortunately, with the present 
data we cannot go much further in this analysis. 

No further improvement is found by adding a third absorp- 
tion region. 

We tried to fit the warm absorber also with "home-made" 
tables built using the CLOUDY photoionization code. Results are 
similar and are therefore not reported here. 

As a consistency check, we then fitted the RGS data with the 
best fit model (all parameters fixed) described above. Spectra 



were rebinned in order to have at least 10 counts per birQ. The 
fit is reasonable (;^^/d.o.f.=333. 4/278), with several, but not indi- 
vidually significant features apparent in the residuals, apart from 
the Ne IX line already included in the spectrum. An improve- 
ment (;\f^/d.o.f. =296. 2/269) is found letting the model parame- 
ters free to vary. The best fit parameters, however, are consistent 
within the errors with those found in the EPIC-pn analysis, and 
the abovementioned residuals are still there. 

We then let the velocities of the absorbing zones free to vary, 
to account for possible inflows/outflows. A better fit is found 
Oir2/d.o.f.=28 1 .2/267), with an inflow velocity of 3900+^^™ km/s 
for the low ionization absorber and an outflow velocity of 
2650^g^g" km/s for the high ionization absorber (see Table [3] 
where "low" and "high" refers to the low and high ionization 
components, respectively, and negative velocities indicate in- 
flows). 

Back to the EPIC-pn analysis, a slight improvement {ls.x^=- 
8) is also found substituting the soft power law with a disk 
thermal component (DISKBB model in XSPEC). The tempera- 
ture is around 160 eV, but the normalization value is unphysi- 
cal, corresponding to about 4x10^ km for the inner radius (the 



' As a check that the rebinning was not altering the fit results, we 
fitted the RGS spectra also with no rebinning and the Cash quality-of- 
the-fit statistics, finding no significantly different results 
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Table 1. Iron line results 



Table 2. Warm absorber results 



Parameters 



Obs. 1 Obs. 2 



Narrow Gaussian 

F„„, (10-"^ cgs) 
EW (eV) 



1 17+0.36 

^•^ '-0.37 

85 

77.2/73 



95+°" 

72 
113.3/93 



Broad Gaussian 






o-(eV) 


143+66 


116!- 


F„„, (10-5 cgs) 


2 1 8+°-^^^ 

^•iO_0.6t) 


1.51!°-^^ 


EW (eV) 


162 


119 


^2/d.o.f. 


66.2/72 


91.4/92 


Diskline 






rin (rg) 


52+'*° 
20!] 


>80 


e 


99+20 


F„„, (10-5 cgs) 


2 24+i-''i 


1 32+°-2' 


EW (eV) 


175 


113 


^Vd.o.f. 


62.6/71 


90.7/91 


Diskline + Narrow Gaussian 






F/,„.,rf,jjr(10-5 cgs) 


2.66!?:°° 


1 4Q+0.57 


EWrf,,, (eV) 


208 


127 


Fli„e,NG (10 5 Cgs) 


0.68!°*? 


0.64!»-^^ 


EW„G (eV) 


44 


52 


^Vd.o.f. 


66.7/72 


96.0/92 



Schwarzschild radius for the 100 million solar masses black hole 
estimated for Mrk 704, Wang et al. 2009, is about one thousand 
times larger). 

A slightly worse fit (A;\f^=14) is instead found if an ionized 
reflection model (REFLIONX, Ross & Fabian 2005) is adopted 
(with or without relativistic blurring) for the soft X-ray compo- 
nent. The power law index is pegged to the lower possible value, 
1.4. 

A moderate improvement (A;t'^=-22) is found letting the ab- 
sorbers be partial. Best fit covering factors of 1 (but loosely con- 
strained) and 0.85 are found for the colder and warmer ionization 
zones, respectively (see Table |2]i. All other parameters are sim- 
ilar to those found in the full absorbers fit. The best fit model 
and data/model ratio can be found in Fig. |5] In Fig. |6] the con- 
tour plots (ionization parameter vs. covering factor) are shown. 
For the low ionization absorber, a full covering is always accept- 
able, while at the 99% confidence level a very low ionization 
(almost neutral) partial absorber is possible. The two parame- 
ters are much better determined for the high ionization absorber, 
with no significant correlation between them. 

Regarding the RGS, a check similar to that described above 
for the full covering scenario was also performed in the par- 
tial covering case. With all parameters fixed, the fit is worse 
than for the full absorbers 0t-2/d.o.f.=356.5/278). A sHghtly 
better fit is instead found letting the parameters free to vary 
(;t'^/d.o.f. =283. 7/267), again with the best fit parameters consis- 
tent within the errors with those found in the EPIC-pn analysis. 



Parameters 


Obs. 1 


Obs. 2 


1 Zone 








10.0"°f' 


4.55+0.42 
L87+}| 


N„ (10^2 cm-2) 


5 16+8-^^ 

3.iO_0 34 


1 5+° °2 


^ (cgs) 


122+' 


39!^ 


E,„„ (keV) 


900+° °°5 


- 


F,,„, (10-5 cgs) 


4 17+0.45 


- 


;tVd.o.f. 


295.1/208 


421.9/230 


2 Zones 






r. 


5.4!|t 


T j:q+0.21 
J.3J_Q14 


r. 


l-64!°;S? 


1-79!°:°^ 


N«,, (10^2 c^-2^ 


0.14!°;»« 


0.085!°™^ 


fi (cgs) 


0.39!°-5^ 


i 1+0.09 


N«.2 (10^2 cm--) 


6.85!°;5^ 


1.63!ii5 


^2 (cgs) 


139!;^ 


1750!?°5 


Efi„. (keV) 


0.904!°;»°5 


- 


F„„, (10-5 cg5) 


A C7+0.49 


- 


^Vd.o.f. 


256.4/206 


279.4/228 


2 Zones, partial covering 






r. 


6 7+^-^<*' 

"•'-3.2 


2 71+0.14 


r* 


1.71+°°^ 


1-18+"^ 


N«., (10^2 cm-2) 


0.10+'™ 
29+^-'^^ 


75+008 
"■'-'-o.oe 


^1 (cgs) 


2.24+0-38 
56+^-^^ 
13.3!^f 
122+jf 
0.38!j|;!J 


Covering Factor (1) 


>0.45 


N;y.2 (10^2 cm-2) 


8.8!1:!1 


ft (cgs) 


97!'5 
0.84!|j;°^ 


Covering Factor (2) 


E,„„ (keV) 


0.895!°°li 




F,„„ (10-5 cgs) 


2 50+°'"* 

^■-'"-0.75 


- 


X^IA.o.i. 


234.3/204 


230.2/226 



(*) The asterisk indicates that a parameter is pegged to one of its 
limits. 



Table 3. Velocities of the warm absorber regions 



Models 



Obs. 1 



Obs. 2 



2 Zones 

V/dw 



-3900!"°° km/s -190±210 km/s 
+2650!^™ km/s -100±250 km/s 



2 Zones, partial covering 

V/ovr 
^ high 



-580±340 
+3000!;50° km/s -^140±300 km/s 
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Fig. 7. RGS best fit spectra and data/model ratio for the first ob- 
servation. Fit with partial covering absorbers. 



The feature-like residuals are still apparent (see Fig.|7]l. Letting 
the velocities of the absorbers free to vary, we found a marginal 
improvement of the quality of the fit (;t^^/d.o.f.=279. 2/265), with 
an outflow velocity of 3000;^[gyQ km/s for the high ionization 
absorber, while the velocity of the low ionization absorber is ba- 
sically undetermined (Table [3]). 

The value of F/, is harder than found in the iron line analy- 
sis. This, together with the relatively high x^, may indicate that 
the warm absorber model we are adopting is not fully adequate 
(e.g., the absorption regions may be more structured that single 
ionization zones, or metal abundances be different than solar). 

For simplicity, in the broad band fitting the iron line is mod- 
eled with a broad Gaussian. No significant changes in the line 
parameters are found with respect to the hard X-ray band fits 
described in the previous section. 

The 0.5-2 (2-10) keV observed flux is 1.53(10.1)xlO-i2 erg 
corresponding to a 0.5-2 (2-10) keV luminosity of 
1.17(2.34)xl0'*-' erg s"', after correction for absorption. 



-2 -1 

cm s 



3.2.2. The second observation 

Due to the much weaker absorption and the absence of the Neon 
emission line, and despite the better statistics, the analysis of 
the second observation resulted to be easier. Also in this case, 
however, the addition of a second zone is definitely required. 

Inspection of the data/model ratio (see FigO shows however 
that some features are still remaining, even if that at about 1.8-2 
ke V is very likely of instrumental origin. Particularly prominent 
is the feature around 1 keV, again suggesting a possible problem 
in the fitting of iron transitions. The addition of a third absorbing 
zone does not cure this problem. 

Also for this observation we find similar results with the 
CLOUDY-based, home-made tables. A significantly worse fit is 
found, instead, substituting the soft power law either with a disk 
thermal component or with an ionized reflection model. 

We then fitted the RGS data with the best fit model (all pa- 
rameters fixed) above described. Given the much better statis- 
tics with respect to the first observation, spectra were rebinned 
to have at least 50 counts per birQ The fit is reasonable 
(;if2/d.o.f.=856.8/679). An improvement 0^^2/(1.0.^=814.4/672) is 
found letting the model parameters free to vary. The best fit pa- 



rameters, however, are consistent within the errors with those 
found in the EPlC-pn analysis. Inspection of the residuals indi- 
cates the possible presence of features, but the only significant 
one (at the 99% confidence level) is an emission line with a cen- 
troid energy of 1.1 14+0.002 keV - to be possibly identified with 
Fe XVll L lines - with a flux of 1.80(+0.75)xlO-5 ph/cm^/s"'. 
However, no such line is found in the EPlC-pn data (upper limit 
of 3.6x10"^ ph/cm^/s '). Finally, no significant inflow/outflow 
velocity is found, with upper limits of a few hundreds km/s 
(TableO. 

A significant improvement (A;t^2=-49) is found letting the ab- 
sorbers be partial. Best fit covering factors of 0.56 (but loosely 
constrained) and 0.38 are found for the colder and warmer ion- 
ization zones, respectively (see Table |2|i. Both absorbers are 
now much thicker and slightly less ionized. The power law in- 
dices are much harder; in particular, the F;, is now 1.18 (even 
if loosely constrained), a value unusually hard for Seyfert galax- 
ies. However, it must be noted that fixing r/,=l .8 a good fit is still 
obtained (;^'2=243. 0/227), with the parameters of the less ionized 
absorber almost unchanged while the more ionized absorber gets 
still more ionized and thicker {Nh,2 ~4x102-' cm"-, ^2 ~230, 
covering factor of 0.27). The best fit model and data/model ratio 
can be found in Fig. |8] In Fig. |9] the contour plots (ionization 
parameter vs. covering factor) are shown. For the low ionization 
absorber the two parameters are almost uncorrected, while for 
the high ionization absorber they are anticorrelated, indicating a 
certain degree of degeneracy in the model. 

The usual check with the RGS is then performed. Fits are 
better than those obtained with the full absorbers, reflecting the 
same improvement found in the EPlC-pn spectrum. With all pa- 
rameters fixed, a ;^'2/d.o.f.=82 1.6/679 is found, while letting the 
parameters free to vary, a;i'2/d.o.f.=752.0/670is found. The best 
fit parameters are consistent within the errors with those found 
in the EPlC-pn analysis. The spectrum and data/model ratio can 
be seen in Fig. [10] Finally, no significant inflow/outflow veloc- 
ity is found for the high ionization absorber (upper limits of 
a few hundreds km/s), while for the low ionization absorber a 
marginally significant inflow velocity of 580+340 km/s is ob- 
tained (Table |3]l. 

For the second observation, the 0.5-2 (2-10) keV observed 
flux is 7.62(1 1.2)xl0"'2 erg cm"^ s"', corresponding to a 0.5-2 
(2-10) keV luminosity of 2.34(2.26)xl0'*^ erg s"', after correc- 
tion for absorption. 

During the revision of the present work, we became aware of 
a paper by Laha et al. (20 11) dealing with the analysis of the sec- 
ond XMM-Newton observation of Mrk 704. The results of their 
analysis are qualitatively similar to ours, the quantitative differ- 
ences probably mostly related to the somewhat different spec- 
tral model adopteqj. The most important qualitative difference 
is that Laha et al. found the two warm absorbers both outflow- 
ing, while in our analysis no significant inflow/outflow is found, 
with only marginal evidence for an inflow of the low ionization 
absorber in the partial covering scenario. 



^ Also in this case, fits to the unbinned spectra provide substantially 
the same results 



^ The main difference in the modeling is that they included a par- 
tial, neutral intrinsic absorber, while keeping full coverage for the ion- 
ized absorbers. If we did the same for the two full absorbers model, a 
slightly better fit is found, which however is significantly worse than the 
one with two partial ionized absorbers. If we did the same for the first 
observation, a fit very similar to that with two partial ionized absorbers 
is found, not surprisingly as the lower ionization zone is almost neutral 
in that observation (see Fig.[6l(. 
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Fig. 8. Obs 2: Best fit model and data/model ratio for the two zones model. Left paneL: full covering absorbers. Right panel: partial 
covering absorbers. 
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Fig. 9. Obs 1 : Contour plots (ionization parameters vs covering factors) for the low (left panel) and high (right panel) ionization 
absorbers. The three curves represent the 67%, 90% and 99% confidence levels. 



3.2.3. Comparison between the two observations 

A comparison between the results of the fits to the two spectra of 
the two observations show, as obvious from Fig. [3] stronger over- 
all absorption in the first observation. Assuming full covering ab- 
sorbing zones, for both of them the column density decreases (by 
a factor about 2 for the colder absorber, a factor about 4 for the 
warmer absorber) and the ionization parameter increases from 
the first to the second observation. If the absorbers are allowed 
to partially cover the source, the colder absorber results thicker, 
warmer and with a lower covering factor in the second observa- 
tion, while for the warmer zone the most significant difference is 
the covering factor, again much lower in the second observation 
(the column density and the ionization parameter both slightly 
increase). 

The detection of a Neon emission line in the first observa- 
tion only cannot be entirely ascribed to the stronger absorption 
(which would make by contrast a constant line more visible). In 
fact, simulations show that a line with that flux should have been 
detected also in the second observation (even if the upper limit 
to such a line, about 10"'', is marginally consistent with the value 
found in the RGS in the first observation). Therefore, an intrinsic 
variability of the line intensity is required, possibly related to a 
different covering factor of the absorbing/reflecting clouds. 



So, the question is: what drove a so large variation of the 
warm absorber over three years? The hard 2-10 X-ray luminos- 
ity stayed approximately the same, while the soft 0.5-2 keV lu- 
minosity was higher by a factor 2 in the second observation 
(even if the reconstructed values, given the large absorption, are 
rather uncertain), which may - in part at least - explain the higher 
ionization parameters. However, the different values in the two 
observations of the column densities and of the Neon line sug- 
gest that also variations in the general properties (size, distance, 
covering fraction, etc.) of the circumnuclear medium have oc- 
cuiTed. 



3.2.4. Time dependent analysis of the second observation 

Due to its intra-observation spectral variability (see Fig. [U, we 
divided the second observation in two parts, corresponding to 
the first 40 ksec (obsA here-in-after) and to the remaining part 
of the observation (obsB), respectively. The three spectra (obsA, 
obsB and total) are shown in Fig[TT] A reasonable fit (y^/d.o.f. 
of 266.2/228 and 274.5/235 for obsA and obsB, respectively) is 
obtained by imposing all spectral parameters to be fixed to those 
found in the total spectrum (two absorbing zones model), leav- 
ing as free parameters only the normalizations of the two power 
laws. The observed spectral variability is then explained mainly 
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Fig. 10. RGS best fit spectra and data/model ratio for the second 
observation. Fit with partial covering absorbers. 

in terms of a variation of the normalization of the soft power law. 
Other choices, like leaving instead as a free parameter the col- 
umn density or the ionization parameter or the covering factor of 
the warm absorbers do not provide acceptable fits. 

Leaving all main parameters free to vary, better fits are of 
course obtained (;t'^/d.o.f=253. 5/221 and 265.2/228 for obsA and 
obsB, respectively) but the parameter which suffered the largest 
variation is again the normalization of the soft power law, con- 
firming the previous analysis. 

4. The Swift and ASCA observations 

The source was observed by Swift five times. In four occasions, 
it was also in the XRT field of view (see Table [Hi, while in the 
fourth observation it was outside. The four XRT spectra, along 
with the two XMM ones, are shown in Fig. [12] 

From the figure, spectral and flux variability is apparent be- 
tween Swift observations (even if far less dramatic than that be- 
tween the two, 3-years apart HMM-Newton observations), im- 
plying a variability time scale of a few months or less. In Table|4] 
the results of the spectral fits are summarized. Due to the modest 
quality of the spectra, a good fit is already obtained with a single 
power law absorbed by a single zone, fully covering absorber 
(no improvement is found letting the covering of the absorber be 
partial). For the lesser quality first and last observations, we fixed 
the power law index to 2, a value similar to those found in the 
second and third observations. The quality of the spectra is in- 
sufficient to establish if the variability is due to a variation of the 
properties of the warm absorber, apart from the first observation, 
when a much larger absorption is apparent. 

Mrk 704 was also observed by ASCA on May 12, 1998, for 
about 35 ks exposure time. A good fit (Y^/d.o.f=348. 9/355) is 
already obtained with a single power law absorbed by a single 
zone, fully covering absorber (see Fig. [T3T l. The source was in 
a lower flux state (2-10 keV flux of about 5x10"'^ erg cm"- 
s"'. The absorber has a column density of 5.6(^j4)xl0^'' cm"^, 
larger than in any other observation, and an ionization parameter 
of 1430(^2™) (cgs). 

5. Discussion and conciusions 

The Seyfert 1 .2 galaxy Mrk 704 was observed twice by XMM- 
Newton, about three years apart (first observation in October 
2005, the second in November 2008). While the hard (>5 keV) 



spectrum was almost constant between the two observations 
(with a possible broadening/strenghtening of the iron line in the 
first observation), the soft part dramatically changed, showing 
much heavier absorption in the first observation. A strong Neon 
IX line is also clearly present in the first observation, while it 
is not detected in the second one. The very similar hard X-ray 
spectrum in the two observations (the behaviour of the intrin- 
sic soft X-ray spectrum is more uncertain because of the strong 
absorption), together with the change in the column density, sug- 
gest that the variations are due to a change in the properties of 
the absorbing clouds. Dramatic changes of the column density 
and/or the ionization state of the absorber have already been ob- 
served in other sources, the most spectacular being NGC 4151 
(e.g. Schurch & Warwick 2002, de Rosa et al. 2007), NGC 7582 
(Bianchi et al. 2009) and NGC 1365 (Risaliti et al. 2005, 2009). 
In the last two sources, in particular, large variations of the col- 
umn density of the (cold) absorber occur on time scales as short 
as less than a day. Here we observe also a strong variation of 
the ionization state of the absorbers. As the primary continuum 
(at least the hard one) stayed almost constant, this indicates a 
variation of the location and/or of the density of the clouds. 

It is interesting to note that Mrk 704 is a "polar-scattering" 
Seyfert 1, i.e. a source with optical polarization aligned perpen- 
dicularly to the radio source axis, as usually found in Seyfert 2s 
(Smith et al., 2004). Smith et al. suggest that in these sources the 
nucleus is seen through the edge of the torus. 

We fitted the spectra with two absorbing regions, either fully 
or partially covering the primary emission. The improvement in 
the quality of the fit when letting the aborbers be partial is only 
moderate for the first observation, but more significant in the 
second one. As the results are somewhat different in the two 
cases, let us discuss them separately. 

With the full covering absorbers, both absorbing zones are 
found more ionized and with a lower column density in the sec- 
ond observation (even if the difference in the column density 
for the low ionization zone is only marginal). The RGS analysis 
of the first observation suggests a possible (but only marginally 
significant) inflow for the low ionization, low column density 
absorber, and an outflow velocity (very poorly determined, with 
90% confidence level range between 1800 and 12000 km/s) for 
the high ionization, high column density absorber. An ionized 
and unstable torus wind, as suggested by Smith et al., may in- 
deed provide an explanation for the latter zone. No significant 
inflow/outflow is instead found, for both absorbers, in the sec- 
ond observation. 

The partial covering scenario, even if providing better fits, 
is more demanding geometrically, requiring a size of the ob- 
scuring clouds of the same order of the size of the emitting re- 
gion. This naturally points to a much closer location of the ob- 
scuring clouds, as e.g. due to a radiatively-driven accretion disk 
wind (e.g. Proga 2003). Interestingly, the variability behaviour 
of Mrk 704 closely resemble that of mini-BAL QSOs (Giustini 
et al. 2010), strengthening the disk wind scenario, at least for the 
more ionized absorber (the colder absorber may be composed 
of orbiting clouds as those found in NGC 1365, Risaliti et al. 
2009). Interestingly, a possible outflow has been detected in the 
first observation (but not in the second, despite very similar col- 
umn densities and ionization parameters). 

The alternative hypothesis, that the partial covering is mim- 
icking the presence of a scattering component originating out- 
side the absorbing region, is ruled out by the short term variabil- 
ity of the X-ray emission. 

Mrk 704 was also observed four times by Swift/XRT be- 
tween January 2006 and January 2007. Spectral and flux vari- 
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Table 4. Swift/XRT results 



Parameters 


Obs. 1 


Obs. 2 


Obs. 3 


Obs 5 


Date 


2006-01-06 


2006-06-14 


2006-09-28 


2007-01-21 


XRT exp. time (s) 


671 


2258 


5602 


1750 


r 


2 (fixed) 


/J_53 


'.Si* 


2 (fixed) 


N„ (10^- cm-2) 
^ (cgs) 


121 

='-'-81 


<1.1 

<450 


Flux (2-10 keV)(10-" cgs) 


0.64 


1.14 


0.88 


0.87 


X^ld.o.f. 


0.31/5 


38.4/36 


49.1/54 


19.3/22 



'% 



iM 



\ 



Energy (keV) 



Fig. 11. The spectra from the first and second part of the second 
'KMM-Newton observation, shown together for comparison. 
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Fig. 12. The spectra of all XMM-Newton and Swift observations 
shown together for comparison. 



ability is apparent between the observations, but the quality of 
the spectra is not sufficient to establish the nature of the varia- 
tions. 

The iron Ka line, studied with XMM, is broad, much broader 
than the optical broad lines but narrower than expected if the 
entire line would be originated in the innermost accretion disk. 
However, if a narrow component (which seems to be almost 
ubiquitous in Seyfert galaxies) is added, then the remaining 
broad component is consistent with emission down to the in- 
nermost stable orbit, even if the quality of the data in not good 
enough to constrain the spin of the black hole. 
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